Nonlinear Fano resonance and bistable wave transmission 
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We consider a discrete model that describes a linear chain of particles coupled to a single-site de- 
fect with instantaneous Kerr nonlinearity. We show that this model can be regarded as a nonlinear 
generalization of the familiar Fano- Anderson model, and it can generate the amplitude depended 
bistable resonant transmission or reflection. We identify these effects as the nonlinear Fano reso- 
nance, and study its properties for continuous waves and pulses. 



I. INTRODUCTION 

The Fano resonance is widely known across many 
different branches of physics; it manifests itself as a 
sharp asymmetric profile of the transmission or ab- 
sorption lines, and it is observed in numerous physi- 
cal systems, including light absorption by atomic sys- 
tems |1| , Aharonov-Bohm interferometer ^ |^] and quan- 
tum dots ^ 1^, 1^ , resonant light propagation through 
photonic-crystal waveguides U 0, M ^| , and 

phonon scattering by time-periodic scattering poten- 
tials 0, IT^ lisl l . From the viewpoint of the fundamental 
physics, the Fano resonance may appear in the systems 
characterized by a certain discrete energy state that in- 
teracts with the continuum spectrum through an inter- 
ference effect. Usually, the discrete state is created by a 
defect that allows one (or several) additional propagation 
paths in the wave scattering which interact constructively 
or destructively. In the transmission line, this interfer- 
ence effect leads to either perfect transmission or perfect 
reflection, producing a sharp asymmetric profile. 

In the classical paper, Fano 0| derived the general for- 
mula, which describes asymmetric line shape of the trans- 
mission or absorption lines. 
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where e = {E—Eji)/{'^q/2) is the dimensionless energy in 
the units of the resonance width 7g, / is the asymmetry 
parameter (Fano factor), and En is the resonance energy. 
In the limit / ^ oo, this formula can also describe the 
so-called Brcit-Wigner resonance profile 0. This quite 
universal formula is used usually to fit a particular pro- 
file observed in experiments and, therefore, to provide a 
proof that the observed phenomenon can be classified as 
the Fano resonance. 

One of the simplest models that can describe the res- 
onant coupling and interaction between a discrete state 
and continuum spectrum is the so-called Fano- Anderson 
model [l3 (see a sketch in Fig. It describes a lin- 
ear " atomic" chain with the nearest-neighbor interaction 
forces and interacting with a defect state through the 



FIG. 1: Schematic of the Fano- Anderson model. An array of 
blue circles correspond to a linear chain and isolated red circle 
is a defect which can be either linear or nonlinear. Arrows 
indicate the coupling between different states. 



nearest neighbors. In application to tight-binding models 
in solids, these forces are linked to a hopping probability. 
This simple model allows to describe the basic physics of 
the Fano resonance in a simple way. In particular, this 
model allows some analytical study that may serve as a 
guideline for the analysis of more complicated physical 
models that predict the Fano resonance effect. 

In this paper we consider an important generalization 
of the Fano- Anderson model under the assumption that 
a single-site defect is nonlinear. We show that this model 
allows to describe a nonlinear Fano resonance, and it can 
generate extremely high contrast between the bistable 
states in its transmission with low input power. We 
study resonant nonlinear transmission and the properties 
of the nonlinear Fano resonance for continuous waves and 
pulses. We would like to mention that a typical physical 
situation when this model can be employed directly for 
describing the resonant effects is the light transmission 
in waveguides and waveguide junctions created in two- 
dimensional photonic crystals with embedded high— Q 
resonators (defects or cavities) with nonlinear response. 
In the nonlinear regime, a defect supports one (or several) 
localized state that is characterized by a discrete eigen- 
value. When the excited localized state interacts with 
the photonic-crystal waveguide, the system can display 
interference effects and Fano resonance. 

The aim of this paper is twofold. First, we study 
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both linear and nonlinear transmission in discrete Fano- 
Anderson models and obtain a number of explicit ana- 
lytical results which allow to get a deeper understanding 
of the physics of this phenomenon. Second, we analyze 
in detail the effect of nonlinearity on Fano resonance and 
demonstrate that, first, it leads to a shift of the Fano res- 
onance frequency and, second, it can lead to bistability 
in the wave transmission. 

The paper is organized as follows. In Section II we 
describe our discrete nonlinear model which provides a 
generalization of the Fano- Anderson model. In Section 
III, we describe the main features of the Fano resonance 
for the case of a local coupling and then study the ef- 
fect of nonlinearity of the defect on the resonant trans- 
mission. In particular, we define the conditions for the 
bistable transmission and show that the nonlinear Fano 
resonance can exist in a broad frequency range. Section 
IV is devoted to the study of a model with nonlocal cou- 
pling in both linear and nonlinear regimes. In Section V 
we consider the wavepacket scattering. Section VI con- 
cludes the paper. 



which allow us to describe elastic scattering process by 
means of a system of nonlinear algebraic equations, 

j 

ojB = EdB + \\B\'^B + VjAj . (4) 

j 

The model Q is the main subject of our analysis of the 
stationary wave scattering. In the following two sections, 
we study first the regime of a local coupling when only 
one coupling coefficient Vo is nonzero, and after that con- 
sider a nonlocal coupling with up to three nonzero cou- 
pling coefficients y_i_o,i- In both these cases, we start 
our analysis from the linear regime and then analyze the 
Fano scattering in the nonlinear regime when A 9^ 0. A 
full time-dependent model (0) is discussed in Section V 
of the paper when we consider the wavepacket scattering. 

III. LOCAL COUPLING 



II. MODEL 

One of the simplest models which describes the physics 
and the main features of the Fano resonance is the so- 
called Fano- Anderson model 0|. In this paper, we use 
a modified version of these model described by the fol- 
lowing Hamiltonian (see Fig. ^ 
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where the asterisk denotes the complex conjugation. This 
model describes the interaction of two subsystems. One 
sub-system is a straight linear chain with the complex 
field amplitude (/)„ at site n which are coupled by nearest- 
neighbor coupling C; this is the system characterized by 
the frequency band of the continuum spectrum. The sec- 
ond subsystem consists of an additional discrete state d 
with the local energy value Ed- For the nonlinear model, 
we assume that the defect possesses a cubic nonlinear 
response, A being the nonlinear parameter. The interac- 
tion between these two subsystems is described by the 
coupling coefficients Vj . 

From the lattice Hamiltonian, we derive a system of 
coupled nonlinear dynamic equations 

i0„ = C(0n-1 + + dy^ VjSnj , 

3 

id = Edd + X\d\^d + YVj(t)j, (2) 

3 

where the dot stands for the derivative in time. For fur- 
ther analysis, we look for stationary solutions of this sys- 
tem in the form. 



First, we analyze the case of a local coupling when 
Vj = Q for J 7^ and Vo is nonzero and study separately 
the linear (A = 0) and nonlinear (A 7^ 0) regimes. 



A. Linear scattering 

We consider the scattering of plane waves with the 
dispersion ujq = 2C cos q propagating along the linear 
chain. Using the second equations of the system Q), we 
find a simple link between two defect-site characteristics 



B = 



VoAo 



LU ~ Ed 

and obtain a single equation 

UjAn = C(An-i + An+l) + 



u - Ed 



(5) 



(6) 



with a one-site scattering potential. For the scattering 
problem, we consider the following boundary conditions 



An — 



/e"3" -Hre-'«", n<0, 
te'^", n > 0, 



(7) 



0„(r) = A„ 



d{T) = Be- 



(3) 



where /, r, and t have the meaning of the incoming, 
reflected and transmitted wave amplitudes far from the 
defect site. Without loss of generality, we will assume 
that incoming amplitude / is real. According to Eq. |(HJ, 
strength of the scattering potential depends on the in- 
coming frequency w, and the system should demonstrate 
resonant scattering. If the frequency of the defect is 
placed in the propagation frequency band of the chain, 
i-e. \Ed\ < 2C, the scattering potential in Eq. ijHJl be- 
comes infinitely large a,t w = tup = Ed, and this will lead 
to total reflection of the incoming wave. 
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FIG. 2; Transmission coeSicient Q calculated for the linear 

Fano- Anderson model ^ in the case Ea — 0, C — 1, and FIG. 3; Frequency spectrum of the system @ vs. the cou- 
A = 0, for several values of the coupling coefficient Vq. pling parameter Vo for the case C = 1, Ed = 0, and A = 0. 



Using the well-known transfer-matrix approach, we ob- 
tain the analytical result for the transmission coefhcient 
defined asT^\t\^/P 



T = 



(8) 



where aq — Cq{ujq — Ed)/VQ and Cq — 2Csing. This re- 
sult corresponds to a simple physics: For the resonant 
frequency ujf, there exist two scattering channels, with 
and without the defect state excited, and destructive in- 
terference between the waves passing these channels leads 
to a complete suppression of the wave transmission (see 
Fig. 121). As a matter of fact, the result (O corresponds 
to the Fano formula , where aq plays a role of the di- 
mensionless energy, Ed is the resonant energy eigenvalue, 
and the resonance width is defined as 



r = 



1^2 



C sinqp ' 



(9) 



where qp is a wavenumber at the resonance, LUp = uj(qp). 
According to this result, the width is proportional to the 
coupling strength V, and it depends on the position of 
the Fano resonance with respect to the boundary of the 
spectrum p^ . In this case, the asymmetry parameter / 
vanishes, and the transmission profile is symmetric. 

Transmission vanishes when aq = 0; this happens at 
the edges of the continuous spectrum band, q = and 
q ^ n, due to the vanishing of group velocities, and at 
the resonant frequency up which, in the case of a local 
coupling, coincides with the discrete eigenvalue of the 
defect mode, LOp = Ed- If the defect is coupled to two 
or more sites, the resonant frequency is renormalized, 
as discussed below. We notice that when the resonant 
frequency is in the middle of the propagation band, the 
transmission profile is symmetric (see Fig. 12), and in the 



case of a local coupling the Fano resonance is observed 
as the resonance reflection. 

Using that |^oP = TI^ and the relation we can 
calculate analytically the amplitude of the defect state, 
that riches the maximum 
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at the Fano resonance frequency up. 

When the coupling coefficient between the defect and 
the chain vanishes, Vo = 0, linear modes of the system Q) 
are described by the continuous spectrum LOq — 2C cos q^ 
while the defect with the frequency oj = Ed is uncou- 
pled. For a finite coupling Vq ^ 0, the defect generates 
two local modes (or bound states) with the frequencies 
bifurcating from the upper and lower edges of the con- 
tinuous spectrum, \ujl\ > 2C, as shown in Fig. The 
appearance of these local modes does not show any char- 
acteristic feature in transmission close to the band edges. 
For the case Ed = 0, these states are located symmetri- 
cally, because the defect frequency is at the middle of the 
spectrum Uq. 

In general, the frequency of the localized state ojl of 
the one-site potential is given by the expression uj^ = 
4:C^ + where ^ is the strength of the potential, and in 
our case, 



ujj = 4C2 



1^4 

^0 



(lol - Edf 



(11) 



Formally, this is the fourth-order equation in ujl, but 
for Ed inside the spectrum cj^, there exist only two real 
solutions which correspond to two bound states. 



B. Nonlinear scattering and bistability 

Now we assume that the isolated defect is nonlinear, 
i.e. A 7^ 0. Using Eqs. and the continuity condition at 
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FIG. 4: Parameter plane {ciq,"iq) defined by different solu- 
tions of Eq. lO- Filled areas correspond to bistability, a 
straight dashed line at the boundary of the bistability region 
corresponds to the Fano resonance. 



the defect site, /+r = i, we obtain the general expression 
for the transmission coefficient. 
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where a; is a real solution of the cubic equation 



{x^ + l){x - aq) - 7q = 0, 



(12) 



(13) 



with the parameter 7^ 

The transmission coefficient defined by Eq. H12() cor- 
responds again to the general Fano formula During 
our derivation, we put x = tan^, where = Arg(r) is 
a phase of the reflection amplitude. Due to the local 
range of the scattering potential it can be shown 
[T^ that the phase of the reflection amplitude is equal 
to 6 — S{q) + tt/2, where 6{q) is a scattering phase shift. 
Therefore, we can consider ea. (|13|) as a nonlinear equa- 
tion for the scattering phase shift, with x — cot (5(g). 
Note here, that ea. H13(l may support complex solutions 
for scattering pha se shift, which will correspond to inelas- 
tic scattering [I^. Such situation is beyond the scope of 
this paper. The expression p2|l is valid only for real so- 
lutions of ea. H13|l . which correspond to elastic scattering 
process. 

Transmission (|12() vanishes at a; = when aq = — 7^ : 



(14) 



or Cq = 0. 



When 7q ~ 0, there exists only one real solution of 
Eq. H13|) . X ^ aq, it leads to the result ||SJ| obtained above. 
In the nonlinear regime (7^ ^ 0), there exist up to three 
real solutions, the output may become a multi-valued 




FIG. 5: Dependence of the transmission coefficient on (a) 
the wave frequency cuq for a fixed intensity I, and (b) input 
intensity / for a fixed frequency uj, for C — A = 1, Vo = 0.8 
and Ed = 0. Bistability regions (|ag| > \/3) are denoted by 
dashed lines. 



function of the input, and this may lead to bistability of 
the wave transmission . 

To define the regions where such a bistability transmis- 
sion can occur, we present the left-hand side of Eq. H13|) 
in the form F(x) = f{x) — jq, so that the bistability 
regions are determined by local extrema of the function 
f{x). When \aq\ < -^/S there is no local extrema, and for 
any value of 7^ (and also A or /) there exists only one 
real solution of Eq. H13() and therefore the transmission 
is always single- valued. For \aq\ > \/3, a pair of mini- 
mum and maximum points of the function f{x) appear 

{aq ± 



at X 



min.max — v"9 \/ "9 ~ 3)/3, respectively, and the 



bistability region is determined by the condition, 

f{Xmin) <lq< .f{Xmax), (15) 

filled on the parameter plane (aq,"fq) in Fig. 01 
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In the nonlinear regime, the transmission coefRcient 
depends on two parameters: the wave frequency ujq and 
input intensity /. First, we study the dependence of the 
transmission coefficient on the wave frequency. From the 
second equation of Eqs. Q), we can see that nonhnearity 
shifts the frequency of the locaHzed mode and, therefore, 
it shifts the position of the Fano resonance. By denoting 
the left-hand side of the Eq. H14|) as g{(^q), it follows that 
g{u}q) is a quadratic function of uiq. If we put Ed — 
then g{—2C)g{2C) < and, therefore, there is an unique 
solution of g{ujq) inside the interval [— 2C, 2C]. In the 
limit of large intensity, i.e. for / — > oo, the solutions of 
the equation g{ujq) = approach the values ±2C. Thus, 
for any value of the input intensity / there always exists 
a single Fano resonance [see Fig. EJa)]. 

To define the position of the Fano resonance analyt- 
ically, we solve Eq. (|14(l explicitly with respect to the 
frequency Wq, 



IV. NONLOCAL COUPLING 

As a matter of fact, the local coupling discussed above 
provides a relatively simple approach for the derivation of 
analytical formulas and understanding the major physics 
of the Fano resonance. However, very often the physi- 
cal problems where the Fano resonance is observed are 
more complicated. In particular, one of the major prop- 
erty of realistic physical models such as photonic-crystal 
waveguides coupled to localized defect modes, is long- 
range interaction and nonlocal coupling 18]. Below, 
we study the generalized discrete (linear and nonlin- 
ear) Fano- Anderson model where the defect is coupled 
to three nearest neighbors of the chain, i.e. in Eq. (0) 
yj^Oforj = -l,0,l. 



A. Linear scattering 



UJp 



2A/2 



(16) 



A proper sign in Eq. H16|) should be chosen depending on 
the particular values of the parameters. 

Thus, by increasing the input intensity, we shift the po- 
sition of Fano resonance such that for / > the trans- 
mission coefficient becomes multi- valued and a bistability 
region appears. Because both aq and 7g depend on the 
frequency Wq, a change of the frequency corresponds to 
a move in the parameter space (aq,7g). For Ed — 0, aq 
changes in the interval [— a™^'^, a™'*'^], where 



(17) 



and bistability is expected for a™""^ > VS [see Fig.EIa)]. 

To study the dependence of the transmission coefhcient 
on the input intensity /, we solve Eq. (|14|l and obtain 



^2 ^ [Uq - Ed)VS 



A(4C2 



(18) 



Equation (|18(l indicates that, depending on the sign of 
nonlinearity A, the Fano resonance appears for any fre- 
quency from the interval [— 2C, or [Ed,2C], and the 
frequency should satisfy the condition. 



In order to find the transmission coefficient for a non- 
local model, we modify the boundary conditions, (O 



An 



/c"?" +rc-'9" , n < -1, 
te'9" , n > 1. 



(20) 



In addition, we should solve a set of coupled equations 
for the sites n = —1,0,1 with the boundary conditions 
(|20|l with respect to the reflection coefficient R = 



Re(5,)2 



Re{bqf + [Cq{u:q-Ed)+lHbq)f 



(21) 



where hq = V^^-f V^ff^Vf + 2e2'«F_iyi + 2e'«Vb(T^-i-hl^i). 
The transmission coefficient T can be found from the 
identity T + i? = 1. 

In the nonlocal case, the transmission coefficient dis- 
plays more complex behavior as compared with the case 
of a one-side coupling (see Sec. Ill above). In particular, 
in this case the Fano resonance is asymmetric. Perfect 
reflection occurs when the second term in the denomi- 
nator of Eq. vanishes. The corresponding condition 
can be rewritten in the form 

Cq{C[LOq - Ed) + + ^i) + 2cosgFiy_i] = 

and it is satisfied at the spectrum band edges (7 = and 
q — TT, due to zero group velocities and also when 



{Ed - LOq)X < 0. 



(19) 



On the plane of parameters {aq,"fq), only jq depends 
on the input intensity / , so that varying / we move along 
the line aq = const and 7^ < or 7^ > (depending 
on the sign of the nonlinearity parameter A). In order 
to achieve bistability, we should satisfy the conditions 
\aq\ > VS and aq^q < 0. The latter condition can be 
reduced to the inequality (Ed — ujq)X < 0, which defines 
the existence of the Fano resonance (|19f) . In other words, 
bistability appears only simultaneously with the presence 
of Fano resonance, see Fig. \^h) and Fig. ^ 



C^Ed~CVo{V-i+Vi) 



(22) 



Equation l|22(l shows that, in the case of nonlocal cou- 
pling to a defect state, the position of the perfect reflec- 
tion gets shifted. If the coupling is strong enough, the 
resonance can move outside the spectrum band. 

In addition, there exists the frequency when a perfect 
transmission occurs. This happens when Re(6g) = or 



C 



'Vo{Vi + V-i)±y/(yf 



4ViF_i)(Fi - V- 



2ViV- 



-<23) 



CO 



FIG. 6: Transmission coeSicient T vs. frequency uoq for the 
three-site interaction, for different values of Vi = Vli. Other 
parameters are: C = 1, Vb = 3 and Ed = 0. According to 
Eq. 1251 . the perfect transmission occurs inside the spectrum 
at Vi = 1.5. There is an interval for coupling coefficient Vi, 
when there is no perfect reflection inside the spectrum. For 
larger values of Vi, both perfect transmission and perfect re- 
flection move to the middle of the spectrum uOq = 0. 

and the corresponding frequency uj\t=i does not depend 
on the discrete state energy Ed- This means when the 
position of perfect reflection lo\t=q is shifted, the point of 
the perfect transmission remains unchanged. This prop- 
erty allows us to vary the width of the asymmetric Fano 
resonance by changing the energy Ed- 

Formally, there exist two solutions for the perfect 
transmission. But in real systems, such as waveguides 
in photonic crystals, some symmetries of the coupling 
coefficients hold, i.e. Vi — V-^i or V^i — Vq and Vi — 0, 
in the cases of three and two nonzero coupling terms, 
correspondingly. For the first case Vi =V-i, we obtain 

w|t=i = -C{Vo/Vi) (24) 

and therefore only one solution for the perfect transmis- 
sion is possible. For small values of Vi, there is no perfect 
transmission, and it takes place at the boundary of the 
spectrum uj\t=i = — 2C, when 

Vi = Vo/2 (25) 

and it exists for any larger value of Vi. When Vi ^ oo, 
the frequency of the perfect transmission moves to the 
middle of the spectrum, uj\t=i (see Fig.EJ. 

In the nonlocal model, the linear spectrum varies sub- 
stantially as a function of the coupling coefficient Vi , see 
Fig. [3 For small values of Vi , two localized states ex- 
ist outside the linear spectrum band ujq (cf. Fig.O. The 
frequency of the upper localized state goes away from the 
spectrum for all Vi . The lower localized state approaches 
the spectrum band but then deviates again for larger Vi . 



FIG. 7: Frequency spectrum for the nonlocal model vs. cou- 
pling coefficient Vi for the same values of parameters as in 
Fig. |S| The insert shows a zoomed area with the double de- 
generacy of the energy levels near the critical value Vcr- 

Levinson's theorem allows us to connect some 
properties of the linear spectrum (FigUJ with transmis- 
sion coefficient (FigEJ. According to this theorem scat- 
tering phase shift at one of the band edge can be repre- 
sented as 

Siq = 0) = niNb + l/2) , (26) 

where Nb is a number of bound state of the system and 
I is a number of quasi-bound states at the band edge. 
Using the relation between transmission coefficient and 
scattering phase shift 19, 21], it can be shown that T{q — 
0) = for even and T{q = 0) = 1 for odd values of I. 
Therefore, we can expect that some quasi-bound states 
should appear at the band edge, when the condition H25() 
is satisfied. The insert in Fig. [7\ shows that exactly for 
this value of the coupling parameter Vi, the spectrum 
band edge becomes doubly degenerate. In this situation 
I = 1. 

For the second case Vq — V-\ and V^i = 0, we obtain 
the following result 

w|t=i = -2C, (27) 

so that the position of perfect transmission does not de- 
pend on the coupling coefficients and it should take place 
at the boundary of the spectrum band. According to 
Levinson's theorem, this may correspond only when some 
quasi-bound state exists exactly at the band edge, I = 1 
in if^ . 

Figure|Slpredicts the perfect transmission at ujg — — 2C 
for all values of Vq, except the point Vo — Vcr when the 
perfect reflection splits off the spectrum band (see 
Fig. ISJ. In this case, for any nonzero value of Vq there 
exists an upper localized state and the quasi-bound state 
at the band edge lu = — 2C (see Fig. O. In addition 



FIG. 8: Transmission coefficient T vs. frequency uoq for tlie 
nonlocal model with two nonzero couplings for different values 
of coupling Vb = V-i. Other parameters are C = 1 and 
Ed = 0. Perfect transmission occurs at uj-i = —2. When Vo = 
Vcr = %/2, the perfect reflection state split off the spectrum 
band and transmission at = —2 vanishes. 



FIG. 9; Frequency spectrum for the nonlocal model with two 
nonzero coupling terms Vb = V-i vs. the coupling coefficient 
V for the same parameters as in Fig. |H] Insert shows a lo- 
calized state a.t u = —2, which does not interact with other 
modes. The second bound state appears at Vb = Vcr when 
the perfect reflection splits off the spectrum band. 



to this, at Vo — V^r another quasi-bound state appears 
at the band edge of the spectrum. In such a situation 
there are two quasi-bound states {I = 2) and according to 
Levinson's theorem H26II transmission coefficient vanishes 
at the bottom band edge. For larger values of coupling 
Vb a second bound state exists in the system. 



B. Nonlinear scattering 

In the nonlinear regime, the reflection coefficient can 
be found in the following analytical form 



R = 



[Re(5,)]2(l+x2) 



(28) 



where a; is a solution of the cubic equation, 



Xl'^c^yl - (x" + l)Re(6,)2(|a,|2Re(6,)2, -f e,y,) = 0, (29) 

where — Re(a^) + Im(a2)a;, ~ Re(a2)a; — Im(a^), 
a, = + e'«(Fo -f e*«Fi), = |a,P[(a;, - Ei)c^ + 
Im(6g)], and the symbols Re and Im stand for the real 
and imaginary parts, respectively. 

We would like to mention again that the general result 
for the reflection coefficient (|28|l is similar to the Fano 
formula where x can be understood as the effective 
dimensionless energy e and the ratio Re(a2)/Im(a2) is re- 
lated to the asymmetry parameter /. This formula H28(l 
includes all other cases discussed above. The condition 
for the perfect transmission Re (6) = coincides with the 
result H23|) , so that the frequency of the perfect transmis- 
sion does not depend on the input intensity /, nonlin- 
earity parameter A, and the energy of the discrete state 



Ed- The position of the perfect reflection does depend 
on all these parameters, and this gives us with a possi- 
bility to vary the width of the asymmetric resonance in 
the nonlocal nonlinear regime. 

Finally, in Fig. ^] we plot the dependence of the trans- 
mission coefficient vs. the input intensity / for different 
values of the frequency. These results are similar to those 
obtained for the model with a local coupling, presented 
in Fig.Efb). 



V. WAVEPACKET SCATTERING 

The previous sections are devoted to the problem of 
the time-independent plane wave scattering. In order 
to verify how these results can be used for describing 
realistic wave scattering, we should analyze the time- 
dependent scattering of wavepackets. In this section, we 
study this problem numerically in the framework of the 
time-dependent model ||2J). 

We consider the propagation of a Gaussian wavepacket 
through a discrete chain with the defect described by the 
system of equations (|2l in the regime of a local coupling, 
when Yq = 'V and V^„^o = 0- As the input wave, we take 
a Gaussian wavepacket of the following form 



0„ (0) = exp 



(n - n^Y 



exp[-i9^(n - no)], (30) 



where is the carrier wavenumber of the wavepacket, 
Iw is its maximum amplitude, a is the spatial width, and 
7io is the initial position. The wavenumber determines 
the velocity of the wavepacket. 
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FIG. 10: Transmission coefficient vs. intensity I for the non- 
local nonlinear model for y_i=Vb = Vi = l,A = l, C=l, 
and Ed ~ 0- Perfect reflection occurs at the edge of the bista- 
bihty domain (see Fig.]!)!. Dashed lines correspond to bistable 
regions. 



For calculating the transmission and reflection coeffi- 
cients, we use the conservation of the norm, and define 
the coefficients as follows 



T = 



n>a 



n<0 



(31) 



where r* corresponds to some time after interaction of 
the wavepacket with a defect. 

In our simulations, we take Ed = and therefore 
qp = 7r/2. First, we compare numerical results for the 
wavepacket scattering with plane- wave analysis in the lin- 
ear regime (see Fig lll|) . When the spectral width of the 
wavepacket 27r/cr is smaller then the width of the reso- 
nance (O, we have observed a good agreement. 

In the nonlinear regime, our numerical results for the 
wavepacket scattering are summarized in Fig. 1121 for the 
dependence of the transmission and reflection coefficients 
on the input wavepacket amplitude 1^, defined in Eq. H3U|) 
for different values of the carrier wavenumber. For the 
case Qw < Qf, Fig. ^1 indicates the existence of the 
Fano resonance characterized by the resonant deep in the 
transmission coefficient and simultaneous resonant peak 
in the reflection coefficient. When > qp, no Fano res- 
onance exists. These results provide a good qualitative 
agreement with the plane- wave analysis presented above. 

In addition to this feature associated with the Fano res- 
onance, there are observed some regions in Fig. 1121 when 
both transmission and reflection coefficient decrease, i.e. 
their sum does not equal to unity. Usually, the dis- 
crete localized mode becomes excited during the inter- 
action of a wavepacket with the defect, and then it re- 
laxes to the ground state. For these particular values of 
the wavepacket amplitude Iw, both defect site and zero- 




1 qp 
wave number, q 

FIG. 11: Transmission coefficient of the wavepacket in lin- 
ear regime A = for different values of the width a. Other 
parameters are C = 1, V = 0.5, Ed = 0. For given parame- 
ters the spectral width of the wavepacket becomes smaller the 
width of the resonance F = 1/4 for a > 20. It leads to a good 
agreement with the analytical result. (Insert) Transmission 
at the position of the Fano resonance qp decays exponentially 
with the spatial width a. 



site particle of the chain get highly excited with long- 
lived oscillations after the interaction (see Fig. I13() . Such 
a behavior is very similar to the scattering by a delta- 
potential with an excited state |l9j | and it corresponds 
to the inelastic scattering process. In our situation the 
localized state is excited due to the nonlinear interaction. 

For the wavenumber q^ — 0.57 such a situation even 
happens near the Fano resonance, i.e. by varying the am- 
plitude of the wavepacket we can observe different out- 
comes for the dynamics of the discrete state (see Fig. ll3ll . 



VI. CONCLUSIONS 

We have studied the properties of wave transmission 
in the presence of a localized nonlinear defect described 
by the generalized discrete Fano-Anderson model. This 
model provides a simple generalization of the well-known 
problem of Fano resonances to the nonlinear case, and it 
can be employed to describe the so-called nonlinear Fano 
resonance effect. We have studied, both analytically and 
numerically, the main features of the resonant transmis- 
sion for continuous waves and pulses, and we have shown 
that the Fano resonance is observed as a specific feature 
in the transmission coefficient as a function of frequency 
and/or intensity. In particular, the presence of nonlinear- 
ity makes the resonance more robust and it can broaden 
substantially the parameter region where the resonant 
transmission is observed, by adjusting either the wave 
frequency or the intensity to satisfy the resonant condi- 
tions. Moreover, we have found the conditions when the 



1| ' I ' I ' V—L 




)i , \ , \ , \ I _^ 

0.5 1 1.5 2 2.5 

Intensity, I (a.u.) 



FIG. 12: Transmission and reflection coeflicients vs. the input 
amplitude 7™ of the wavepacket defined in Eq. H^O^ for C = 1, 
V = 0.5, Ed = Q, \ = l, a = A, and no = -10. 

Fano resonance is associated with the bistable transmis- 
sion. 

We would like to emphasize that many of the effects 
described in this paper can find applications in a vari- 
ety of different nonlinear physical systems, including the 
transmission of straight and curved waveguides coupled 
to localized defect modes in photonic crystals. Moreover, 
a number of interesting effects already observed in some 
of the nonlinear systems can be identified as the mani- 
festation of the nonlinear Fano resonances. In particular, 
we believe that the resent observation of excitonic opti- 
cal bistability in CU2O (22j can be explained in terms of 




Time 

FIG. 13: Evolution of the defect site during and after the 
wavepacket scattering, for = 0.57. Other parameters are 
the same as in Fig. 1121 Different final states of the defect are 
observed for a slight variation of the amplitude of the incom- 
ing wavepacket, being a signature of the system bistability. 
Inset shows a log-plot of this dependence and indicates that 
the defect state relaxes to the ground state algebraically. 

the nonlinear Fano resonance which could shed additional 
light on the physics of these experimental observations. 
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